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We study the radiative corrections to all K(a decay modes to leading non-trivial 
order in the chiral effective field theory, working with a fully inclusive prescription 
on real photon emission. We present new results for K^ modes and update previous 
results on K e ^ modes. Our analysis provides important theoretical input for the 
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extraction of the CKM element V us from K^ decays. 
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1 Introduction 

With the advent of precision measurements in Kaon physics (see [1] and references therein), 
Kt 3 decays offer the opportunity to probe charged current weak interactions at unprece- 
dented levels. Most notably, with current experimental uncertainties, K^ decays allow 
one to access the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing angle V us at the sub- 
percent level, and also provide competitive probes of lepton universality and the ratios 
of light quark masses. In order to fully exploit the amazing experimental achievements, 
it becomes mandatory to have theoretical control of these decays at the percent level or 
better. This requires quantitative understanding of the vector and scalar K — ► n form 
factors as well as the electromagnetic (EM) corrections. The framework to analyze the EM 
corrections is provided by Chiral Perturbation Theory (ChPT) [2], the low energy effective 
field theory (EFT) of QCD, extended to include the photon [3] and the light leptons [4] as 
active degrees of freedom. ChPT exploits the special role of ir, K, 77 as Goldstone modes as- 
sociated with the spontaneous breaking of chiral SU(3)l x SU(3)r symmetry, and provides 
a systematic expansion of the amplitudes in powers of the masses of pseudoscalar mesons 
and charged leptons (p ~ M n ^ t e/A x with A x ~ 47T.F 7r ~ 1.2 GeV) and the electromagnetic 
coupling (e). 

In this article we present results on the electromagnetic corrections to the four K^ 
decay modes (K = K ± , K°; £ — e, /j,), based on a calculation of the amplitudes to leading 
non-trivial order in ChPT (0(e 2 p 2 )). For all modes, we focus on (i) the electromagnetic 
(EM) corrections to the Dalitz plot, which are needed to extract the momentum dependence 
of the K — > 7r vector and scalar form factors from the experimental distribution; (ii) the 
integrated radiative correction to the decay rate, which is a crucial input in extracting 
the CKM mixing angle V us from K — > 7r£z/[7] decays. The theoretical framework for 
the calculation of electromanetic contributions to 0(e 2 p 2 ) in K^ decays was presented in 
Ref. [5] and full numerical results on the K e3 modes were given in Refs. [6], adopting a 
specific prescription for treating real photon emission and a specific factorization scheme 
for soft photons , which results in the partial inclusion of higher order terms in the chiral 
expansion. The novel features of the present work can be summarized as follows: 

• Rather than using the soft-photon factorization procedure of Ref. [5], we work here 
to fixed chiral order e 2 p 2 , providing the complete corrections to decay distributions 
and total decay rates to 0(e 2 p°). 

• We give new results for K^ modes and update our previous analysis of K e3 modes. 

• We use a fully inclusive prescription for real photon emission, which is more appro- 
priate for comparison with the experimental results. 

• We update the structure-dependent EM correction, using the recent estimates of the 
relevant low-energy constants (LECs) provided in Refs. [7, 8]. 

Preliminary results of our analysis have been made public in conference talks and 
proceedings [9, 10] and should be considered obsolete after the current publication. The 



paper is organized as follows: in Section 2 we give an overview of various contributions 
to Ki^ radiative corrections and derive the relevant master formula for the corrections to 
fixed chiral order (0(e 2 p°)). In Section 3 we present our results for differential and total 
radiative corrections, discussing their uncertainty. In Section 4 we present our conclusions. 

2 Radiative Corrections to K^ decays: overview 

2.1 Generalities on K^ decays 

Let us briefly recall the main features of Kf$ decays. The invariant amplitude for the 
process K(j>k) — ► ^(Pn) ^ + (pe) ^e(Pu) reads 



M = ^C u(p v ) Y (1 - 75) v(p e ) C K 



irit) ( PK + Pw ), + f K -(t) ( PK - ^)m 



(i) 



where Ck = 1 for K® 3 and Ck = l/\/2 for K^ modes. The expression in square brackets 
corresponds to the matrix element (7r(p n )\V^~ l5 \K(pK)) , expressed in terms of the form 
factors f±(t), which depend on the single variable t = (px — P-w) 2 and are known to 
0(j» 6 ) in ChPT [11, 12, 13]. To this order, a number of unknown LECs appear, of which 
only a few can be determined experimentally. A complete prediction to 0(p 6 ) requires 
theoretical input beyond ChPT, either from analytic approaches [14] or lattice QCD [15]. 
For phenomenological applications, it is common to parameterize the form factors /+(£) 
and fo(t) = f+(t) +t/(Mf C — M^)/_(i) in terms of slope and curvature parameters x which 
can then be measured: 

f *' (t) ~WW) = ' + A+ SIT + 2 A + mu • (2) 

/5 . (t ,,£w %^( Ao _ A+ _^< ,. (3 ) 



The spin-averaged decay distribution depends on two independent kinematical vari- 
ables, which we choose to be 

_ 2p w ■ p K _ 2^ _ 2p K ■ p e _ 2E(_ 

Z ~ Ml ~ M K ' V ~ Ml ~ M K ' [ ] 

where E n (Eg) is the pion (charged lepton) energy in the kaon rest frame, and Mr indicates 
the mass of the decaying kaon. Then the distribution (without radiative corrections) reads 

jp(o) r 2 \v I 2 /if 5 r 2 

g_ = ^1^^ ^.(0)1'^,,), (5) 

p""(y,z) = Af(y,z) \f?'(t)\* + Af(y,z) /**(«)/-*(*) + Af(y,z) |/5*(f)| 2 ,(6) 



'"See Ref [16] for a dispersive parameterization of the scalar form factor. 



where the kinematical densities read (re = (me/Mx) 2 and r n = (vti^jMk) 1 )'■ 

Af\y,z) = A(z + y - 1)(1 - y) + r e (Ay + 3z - 3) - Ar n + r £ (r n - r e ) , 

Af(y,z) = 2r e (3-2y-z + r e -r n ) , 

Af(y,z) = r e (l + r w -z-r e ) . (7) 

In the analysis of K e s decays, the terms proportional to A 2 3 can be neglected, being 
proportional to r e ~ 10~ 6 . Finally, the decay rate reads 

n°\K a ) = 9MY^^k | /f ir (0) |2 7 (o) (At) (8) 

I { k\(\) = JdydzfP\y,z), (9) 

where the integral extends on the physical domain T>^ defining the three-body Dalitz plot 
(see Ref. [5] for the explicit definition). 

2.2 Radiative corrections: soft factorization vs fixed chiral order 

The above description of differential distributions and decay rates is modified by EM effects, 
which involve the emission of both virtual and real photons. Short distance electroweak 
corrections can be lumped in an overall factor S ew — 1 + — (l — f 1 -) x log^ + 0(^% a ), 
which is common to all semileptonic charged- current processes [17]. Long distance EM 
corrections to Kg^ decays can be studied within ChPT. The leading non-trivial corrections 
to the amplitudes appear to 0(e 2 p 2 ) and imply corrections to the form factors, differential 
distributions, and decay rates starting to 0(e 2 p°) 2 . 

In Ref. [5], it was argued that long distance EM effects can be taken into account by (i) 
a universal (i.e. non structure-dependent) shift in the densities A\ (y,z) accompanied by 
(ii) structure-dependent corrections to the form factors f± n (t). This result was obtained by 
factorizing out of the amplitude the universal soft photon corrections [18] that are sensitive 
only to charges, masses, and momenta of the particles involved in the decay. While this 
recipe has the benefit of being simple and elegant, it inherently mixes different orders in the 
chiral power counting (e.g. the soft-photon corrections proportional to f^ w (t) only appear 
in the EFT calculation to 0(e 2 p 4 )). As a consequence, the resulting decay distribution 
and rate contain not only the full chiral corrections of order e 2 p° but also incomplete 
higher order corrections, generated by the factorization procedure. Since we are studying 
a fully photon-inclusive rate, there are no large logarithms associated with the factorized 
soft-photon corrections and therefore the partial higher order corrections that are included 
with the recipe of Ref. [5] are not expected to give the dominant contributions to any given 



2 Sincc u(p v ) 7 M (1 — 75) v(pe) ■ (pk ±Ptt)^ ~ 0(p 2 ), EM corrections to f±* start at 0(e 2 p°). Moreover, 
since y, z^rg^r^ ~ O(l) we can book the densities AJ23 as quantities of 0{\). Therefore, corrections of 
O(e 2 p ) to f± v induce corrections to the decay distributions and rates of O(e 2 p ) (see Eqs. 6 and 8). 



order: cancellations with unknown terms are possible. Motivated by this, in the present 
work we give the corrections to decay distributions and rates to fixed chiral order, namely 
0(e 2 p°), to which the complete answer is known. We shall then use the comparison with 
the procedure of Ref. [5] as a validation of our estimate of the theoretical uncertainty. 

The virtual photon corrections to all K^ amplitudes (K = K ± ,K° and £ — e, /i) are 
known to 0(e 2 p 2 ) [5, 6], while the real photon emission was worked out explicitly in those 
references only for K e3 modes (and only for a specific prescription on the treatment of 
real photon emission [19]). Our goal here is to provide a unified discussion of radiative 
corrections of 0(e 2 p°) to all K^ decay rates, working with the fully inclusive prescription 
on real photon emission. We now sketch the derivation of the corrections induced by 
virtual and real photon emission to fixed chiral order [0(e 2 p°)}, and how they combine into 
a master formula for the inclusive rate. 

2.2.1 Virtual photons 

One-loop amplitudes involving virtual photons, together with the associated local coun- 
terterm contributions, induce an effective correction of 0(e 2 p°) to the QCD form factors 
f± n , of the form: 



a 

4:71 



f?*(t) - f^(t) + 5f^(v) + -T c (v,mlM 2 ;M 2 ), (10) 



f?*(t) - f™(t) + 8f?*(v), (11) 

where M c is the relevant charged meson mass and v — u = (px ~Pe) 2 for K^ 1 decays while 
v — s = (Pn+Pe) 2 for K° decays. The function T c (v , raj, M 2 ; M 2 ) [5] encodes the universal 
soft photon virtual corrections and is infrared divergent (thus it depends explicitly on 
the infrared regulator M 7 ). On the other hand, the corrections 5f± n (v ) encode structure 
dependent effects through one- loop corrections and chiral low-energy constants [5]. We 
report their expressions in Appendix A in terms of functions defined in Ref. [5]. 

Keeping in mind the chiral properties of f± n (t), namely that /+(£) = 1 + 0(p 2 ) and 
f-(t) = 0(p 2 ), the effect of virtual corrections to leading order in ChPT amounts to the 
following 0(e 2 p°) shift to the differential distribution of Eq. 6: 

5p EM - vht ™\y,z) = A?\y,z).[2 5f^(v) + ^r c (v,m 2 ,M 2 ;M 2 : 

+ A?\y,z)-8f««(v) . (12) 

2.2.2 Real photons 

It is well known that only the inclusive sum of K — > tt£i> and K — > iriu + w-f (with 
n — 1, 2, ...) decay rates is infrared (IR) finite and observable. In the chiral power counting, 
the leading contribution to the radiative amplitudes K{j>k) —> 7r(p ff ) £ + (pe) ^e.{Pv)l{k) is 
of 0(ep): this is all we need for the analysis of the rates to 0(e 2 p°). To this order the 



radiative amplitudes read: 
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The resulting correction to the differential (or total) decay rate can be calculated from 

Mr, ^ I 



dT{K -> tt£z/ 7 ) 



2M. 



K 



pol 



1 2 dVt n i vl 



d 3 pi 



II iri S (A) (p K - P* - Pe - Pv - k) , 

Pi 



(15) 
(16) 



j=7r,<;,^,7 



where dQ n i uy is the 4-body Lorentz invariant phase space. 

The square modulus of the radiative amplitude can be decomposed into the sum of an 
IR singular term (TjpJ and an IR finite inner bremsstrahlung term (Tib), as follows 
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with the IR singular pieces given by: 
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(20) 



The terms T^/ generate an infrared divergence when integrated over the soft photon region 
of the four-body phase space, while the remaining terms are IR finite. Integrating over all 
variables except y and z, leads to a correction to the Dalitz plot density of the form 



Sp 



EM-real 



l(0). 



a 



(y,z) = A?>{y,z) ■ -I (y,z; M 7 ) + A?(y,z) , 

7T 



(21) 



where Io(y,z; M 7 ) arises from Tir while A* B from Tib- The functions Io(y,z;M 7 ) and 
A^ B (y,z) depend on the cut on hard real photon emission. Results based on integrating 
over all kinematically allowed photon energies are known analytically [5, 6, 19]. 

Finally although both Io(y, z; M 7 ) and T c (v, mf, M%\ M 2 ) are individually infrared di- 
vergent, they combine into the IR finite function 



A m (i/,z) 



a 



I (y,z-M 1 ) + -Y c (v,mlM 2 -M, 



f22) 



and we end up with a finite correction to the Dalitz plot density: 
6p EM (y,z)=A?\y,z)-[A m (y,z) + 2 8f«*(vj\ + A? (y,z) +A?\y, z) ■ <5/5» . (23) 

2.3 Master Formula 

After inclusion of both long-distance and short distance (S cw ) radiative corrections, the 
differential decay distribution reads: 



dY 



(JTf \V n .J Mrs Ciz 



-djdz T 128 7T3 5 - l/^(0)| 2 [p<°>(y,*) + ^ EM (y,^) 



(24) 



This expression should be the basis to properly determine experimentally the momentum- 
dependence of the QCD form factors f±(t) appearing in p^(y, z). 

Corrections to the decay rate are obtained by integrating over the variables y and 
z. In the case of fully inclusive prescription on the radiated photon, the real photon EM 
correction should be integrated not only over the 3-body region T>^ but on the whole region 
£>4 allowed by 4-body kinematics. Taking into account all these corrections, the master 
formula for K^ decay rates reads: 



F (^3[ 7 ]y 



G\\V us \ 2 M\Cl 



(o), 



s e w \fr~m 2 /$(* 



128 7T 3 

where the strong isospin breaking correction is 
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Jf^(0), 
and the EM radiative correction 

<*EM = ^EMps) + <*EMp4-3) 

receives contributions from the 3-body and 4-body kinematical regions: 
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4-3J 
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I K e(^i 



2vr 4 M| JC4 _ 3 



dQireu-y (T m + T IB ) 



(25) 

(26) 

(27) 

(28) 
(29) 



Note that Ginsberg's prescription [19] for real photon emission, which was adopted in 
Refs [5, 6], amounts to discarding the integral in Eq. 29. 




a 



0.05 


-0.05 









Kj 3 






1 

0.9 


\ 


\ 






1 _ 

/■ 


0.8 




\ 




1" 


/ - 


0.7 






x 


w 


- 


0.6 






\l 


P' 


- 



0.1 0.2 0.3 0.4 0.5 0.£ 



0.7 0.8 0.9 



0.2 0.3 0.4 0.5 



Figure 1: Density plot of the EM correction to the differential distribution 
(5p EM (y,z)/pW(\i)(y,z)) of K% (left panel) and K± (right panel). 



3 Results and discussion 

The main outcome of our analysis are the differential corrections 5p EM (y,z) to the Dalitz 
plot density (Eq. 23) and the integrated corrections 5^{V 3 ) and S^iV^^ (Eqs. 28 and 
29 respectively). 5p EM (y, z) is known analytically through the work of Refs. [5, 6, 19]. The 
integration needed to calculate S§^(T> 3 ) has been performed with the Gauss quadrature 
method. On the other hand, the integration needed to calculate S^(V^ 3 ) has been 
performed with a Monte Carlo technique based on the RAMBOS event generator [20]. 

In order to give numerical results, we have to specify the input parameters. The 
differential and integrated EM corrections to 0(e 2 p°) depend on a number of LECs of 
ChPT. The EM LECs are given by convolutions of appropriate QCD correlators with 
known kernels. They have been recently estimated in Refs. [7, 8] by replacing the QCD 
correlators with meromorphic approximants, in the spirit of the large- Nq expansion. We 
use the results of Refs. [7, 8] for our central values, and conservatively assign 100% fractional 
uncertainty to the LECs. 

The integrated corrections S^(T> 3 ) and ^^(^4-3) also depend, through the normal- 
ization factor I K \{\i), on the slope and curvature of the scalar and vector form factors, 
namely A+, A + , and Ao- For these quantities we use the results of a global fit to all con- 
sistent experimental data (i.e. without including the NA48 result) as reported in Ref. [1]: 
A+ = (25.0 ± 0.8) ■ 10~ 3 , A+ = (1.6 ± 0.4) ■ lO" 3 , A = (16.0 ± 0.8) • 10~ 3 . The choice 
of this reference set of input parameters is very simple but somewhat inconsistent as the 
charged and neutral K parameters are distinguished by strong isospin breaking and EM 
effects [13, 23]. However, one should keep in mind that the uncertainty on I K \{\j) induced 
by the Aj is below the percent level and is completely negligible in the analysis of 5^ 
(although it has some impact on the extraction of V us through Eq. 25). In Table 1 we 
provide the I Ki (\i) corresponding to our choice A,, so the reader can easily convert our 
results for 5^ to any choice of slope parameters with a simple re-scaling. 




Figure 2: Density plot of the EM correction to the differential distribution 
(5p EM (y,z)/pW(\i)(y,z)) of K° 3 (left panel) and K% (right panel). 



3.1 Corrections to the Dalitz plot 



In Figures 1 and 2 we show a density plot of the ratio 8p 



^EM/^ 



/p(°)(\) for K%,K%K% and 



K^ as a function of the variables y, z, corresponding to the input on Aj and LECs specified 
in the previous subsection 3 . The theoretical uncertainty on the LECs and higher order 
corrections induces an overall uncertainty of about ±0.3% in Sp EM /p^(Xi). It is important 
to notice that the correction to the Dalitz distribution can be locally large (0(10%)) and 
does not have definite sign, implying possible cancellations in the integrated total EM 
correction. 



3.2 Corrections to the decay rates 

Table 1 summarizes the numerical results of the long-distance radiative corrections to fixed 
order e 2 p°, obtained using the central values for the LECs, slopes, and curvature of the 
form factors as described above. 

Two prominent features of the results in Table 1 can be understood on a qualitative 
level. First, the EM corrections to the neutral K decays are expected to be positive and 
sizable on account of the Coulomb final state interaction term between £ + and 7r~, that 
produces a correction factor of na/vlf ~ 2% over most of the Dalitz plot. While the 
exact correction and the relative size of K® 3 and K® 3 depend on other effects such as the 
emission of real photons, the qualitative expectation based on Coulomb interaction is con- 
firmed by the detailed calculation. Second, the large hierarchy S^iV^s) <C 5^ (^4-3) 
admits a simple interpretation in terms of bremsstrahlung off the charged lepton in the 
final state. The probability of emitting soft photons is a function of the lepton velocity 
V£ which becomes logarithmically singular as ve — > 1, thus enhancing the electron emis- 
sion. For typical values of ve in T>4_3, the semiclassical emission probability [22] implies 

OZ>4-8)/<©(2>4-8) ~ 20 - 40. 

The theoretical uncertainty to be assigned to 5§^ arises from two sources: the input 
parameters (LECs, A*) appearing in the 0(e 2 p°) correction and unknown higher order 



3 Numerical tables for these corrections are available from the authors upon request. 



terms in the EFT expansion, starting at 0(e 2 p 2 ). For the parametric uncertainty we find 
that: 



Experimental errors on the form factor parameters A+, A + , and Ao induce a fractional 
uncertainty in I K \{\) and in 8^ well below the percent level [1]. We can safely ignore 

Ki 

EM- 



this source of uncertainty in S^ 1 



To the order we work, the electromagnetic LECs contribute a ^-independent term to 
8f± n (v), which thus affects the decay rates as follows, 



(0)/„. „\ / j„.j^/|(0) 
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5™(V 3 ) ~ 25f + _ • "* ^— + 6f. _ • ^ ^— . (30) 
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r(0) 



[("'(l.) LECs 7" W CA- 
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r(0), 



The coefficient of 2#/ + is 0(1) while the coefficient of 5/_ is roughly 0.2 for K^ 3 
decays and completely negligible for K e3 decays, being 0(m e /Mx) 2 . Taking a very 
conservative attitude, we assign a 100% fractional uncertainty to 5,/±|lecs (using 
LECs central values from Refs. [7, 8]), which induces an absolute uncertainty in all 
the 6£ M of about ±0.1%. 

In order to discuss the error coming from higher order chiral corrections not included 
in our analysis, we find it useful to decompose (to each chiral order) the EM corrections 
S Ke in terms of Si 2 3 4: 



s-K°e 


= 5! + 5 2 + 5 3 + 5 4 , 


§K°n = 


= 5i + 5 2 - S 3 - 84 , 


5 K±e -- 


= Si - 5 2 + S 3 - 5 4 , 


6 K± » -- 


= Si — S 2 — S 3 + S4 . 



(31) 

Here Si represents a correction common to all modes, S 2 a correction anti-correlated in 
kaon isospin but blind to lepton flavor, and finally S%^ are lepton-universality breaking 
terms, correlated and anti-correlated in kaon isospin, respectively. To 0(e 2 p°) we find 
SfP° = 0.63%, Sf p ° = 0.57%, Sf p ° = -0.08%, S( p ° = -0.12%. On the basis of chiral 
power counting we expect the higher order corrections to scale as 

Sf p2 ~ (M K /(4nF n )) 2 ■ Sf p ° ~ 0.2 • Sf p ° . (32) 

This estimate is validated by comparison of the fixed chiral order results with the ones 
obtained within the "soft-photon factorization" approach discussed in Section 4 of Ref. [5] , 
which include a class of higher order chiral corrections (see Table 2) 4 . The only anomaly 



4 Note that in Section 5.3 of Ref. [5] an alternative factorization prescription is given, which is valid 
only for K e z modes. The latter prescription was used in the numerical analysis of Refs. [5, 6] and would 
lead to results slightly different from those of Table 2. The first two entries in the first column should be 
replaced as follows: 0.41 — » 0.56 and -0.564 — » -0.41. 
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IrA^) 


S(»3)(%) 


^(^4-3) (%) 


0%) 




0.103070 
0.105972 
0.068467 
0.070324 


0.50 
-0.35 

1.38 
0.007 


0.49 
0.45 
0.02 
0.009 


0.99 ± 0.30 

0.10 ± 0.30 

1.40 ± 0.30 

0.016 ± 0.30 



Table 1: Summary of phase space integrals and EM corrections to the K^ decay rates. 
The EM corrections are calculated to fixed order in ChPT (0(e 2 p°)). The phase space 
integrals are calculated using slope and curvature parameters from the fit of Ref. [1]. The 
uncertainty estimate is discussed in the text. 





S^(Vs)(%) 


<^p4- 3 )(%) 


*&(%) 


K% 


0.41 
-0.564 

1.57 
-0.006 


0.59 

0.528 

0.04 

0.011 


1.0 
-0.04 

1.61 
0.005 



Table 2: Summary of EM corrections to the Kjn, decay rates calculated according to the 
"soft-photon factorization" approach of Ref. [5], which includes incomplete higher order 
terms in the chiral expansion. Comparison with the results of Table 1 validates our estimate 
of the theoretical uncertainties. 



appears to be in the coefficient £3, where we find £3 : —0.08% — > —0.16% when going from 
fixed chiral order to the soft factorization scheme. This can be traced back to the cancel- 
lation between the negative contribution from T> 3 (-0.31%) and the positive contribution 
from 224-3 (0.23%). Multiplying these individual pieces by 0.2 gives ~ 0.06 and ~ 0.05, 
respectively, which is just the order of magnitude of the shift we are seeing (-0.08 — ► -0.16). 
Based on the above discussion, we bound the higher order uncertainties as follows: 
\Sf p '\ < 0.13%, \5f p2 \ < 0.11%, \5f p2 \ < 0.08%, \8f p2 \ < 0.025%. Adding these linearly 
we estimate the uncertainties quoted in Table 1 for the total corrections. Finally, using the 
same bounds on 6% p we estimate the theoretical uncertainties on the linear combinations 
which are relevant for lepton universality and strong isospin-breaking tests: 



rK u e r 

Orm - 



'EM 

°EM 

K K±e 
°EM 

°EM 



K°ix 

EM 



°EM 

?K°e 
°EM 

°EM 



-(0.41 ±0.20)% , 
(0.08 ± 0.20)% , 
-(0.89 ±0.30)% , 
-(1.38 ±0.30)% . 



(33) 
(34) 
(35) 
(36) 
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4 Conclusions 

In this work we have provided a unified discussion of the radiative corrections of 0(e 2 p°) to 
all Kiz decay rates. We have argued that through the calculation of K — > -kIv amplitudes 
to 0(e 2 p 2 ) in the chiral effective theory we can derive the complete corrections of 0(e 2 p°) 
to the Dalitz plot density and the integrated decay rate. We have systematically discarded 
higher order effects that are only partially known, and we have included the unknown 
effects in a generous estimate of the theoretical uncertainty. For the first time we have 
presented complete numerical results for the K^ modes, while also updating the previous 
analysis of K e3 modes. 

The main outcome of our investigation is summarized in Table 1, which contains the 
corrections to the total (fully photon inclusive) decay rates of all K& decay modes (K = 
K ± , K°; £ = e, //). These results provide important theoretical input for the determination 
of the product f+ n (0) ■ V us from K& decays at the 0.2% level [1] through Eq. 25. This 
task requires as additional theoretical input the factor 8g-J, 2 x, which has been recently 

updated in Refs. [13, 23]. We refrain here from producing a number for f+ n (0) • V us : 
the result reflecting most recent experimental data can be found in the Flavianet Kaon 
Working Group web page [24]. 

Further reduction of the theoretical uncertainty on the corrections <5™ would require 
an analysis of the amplitudes to 0(e 2 p 4 ) in ChPT, which is beyond the scope of this work. 
At the moment there appears to be no immediate need for such an analysis, since the 
error on V us is dominated by the ~ 1% theoretical uncertainty in f? n (0), for which a 
compilation and discussion of theoretical results can be found in Ref. [1]. 
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A Corrections induced by virtual photons: Sf^^^v 

Using the notation of Refs. [5, 6], and identifying y™~ loc f R e f. [5] with f± of Ref. [6], we 
have: 






a 

47T 

a 

4n 



T^mlMl) + T 2 (u,mlM 2 K ) 
T^mlMl) -T 2 (u,m 2 ,M 2 K ) 



T J+ 1 



, ?K+7r 



(37) 
(38) 
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and 



*/+ 



5f«°«- (s) = ± 



Y*(s,m\,M*)-Y x {s,m\,Mi 



+ lf*~ , (39) 



+ f?°*~ ■ (40) 
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